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Abstract 



Low-energy muon beams are useful for a range of physics experiments. We consider the production of low-energy (i + beams with 
small energy spreads using frictional cooling. As the input beam, we take a surface muon source such as that at the Paul Scherrer 
£S| Institute. Simulations show that the efficiency of low energy \n + production can potentially be raised to 1%, which is significantly 
£H higher than that of current schemes. 
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1. Introduction 
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Intense low-energy muon (LE-U+) sources with small energy 
^ spreads are important for many experiments including mea- 
surements of rare muon decays, fundamental constants, and 
^ the muon's anomalous magnetic moment; determination of the 
Lorentz structure of the weak interaction; and tests of CPT 
' 5d s y mmetr y ar >d QED. They are also useful in condensed matter 
£>~j physics: muon beams with small energy spreads and mean en- 
ergies tunable from eV to keV can be used as sensitive probes to 
, study various phenomena in materials 01 . A muon collider J2l 
and neutrino factory [3] also require high-quality muon beams. 

Muon beams are tertiary beams, produced in the decay of 
pion beams that are produced by proton beams hitting targets. 
They initially have large energy spreads that must be reduced. 
This process, called "cooling," is important for all the exper- 
£f) imental endeavors described above. Frictional cooling is one 
promising method to reduce the energy spread. 

In this paper, we explain the frictional cooling concept and 
detail a scheme for the production of a LE-JJ+ beam from a sur- 
face muon source, such as the one available at the Paul Scherrer 
Institute (PSI), Switzerland fl]. We then compare the efficiency 
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of our cooling scheme to the current LE-|i + 

r$ ]psi & 



2. Frictional Cooling 

Frictional cooling is the bringing of charged particles to an 
equilibrium energy by the balancing of energy loss to a mate- 
rial with energy gain from an electric field [6]. Figure Q] shows 
the stopping power, - j^, as a function of kinetic energy T for 

\i + , where p is the density of the medium and is the energy 




'Corresponding author 

Email address: baoyuSihep .ac.cn (Yu Bao) 



10 5 10° 
T [keV] 



Figure 1: Stopping power for u, + in helium as a function of kinetic energy T. 
T c q is the equilibrium energy. The red line represents the accelerating power of 
an external electric field. The shadowed band is the frictional cooling region. 



loss per unit path length. The stopping power has been calcu- 
lated using velocity-scaled proton data [7]. Applying an elec- 
tric field to restore kinetic energy in the longitudinal direction 
brings the particles to the equilibrium energy (r eq ). Particles 
with kinetic energies less than T eq accelerate because they gain 
more energy from the electric field than they lose to the ma- 
terial. Particles with kinetic energies greater than T eq but less 
than T (see Fig. [TJ decelerate because they lose more energy 
than they gain. Thus all particles with kinetic energy less than 
T' ecj approach T eq and the energy spread decreases. 

Frictional cooling requires particles have low kinetic ener- 
gies. Near T eq , where particles lose energy through excitations 
of the material atoms, charge exchange processes, and nuclear 
recoil, the stopping power of the retarding medium is large. The 
average density must be kept low so that 4p can be balanced by 
realistic electric fields. 

One way to achieve this is using a series of thin foils with 



Preprint submitted to Elsevier 



January 18, 2010 



an electric field present between them, as was used by a previ- 
ous experiment at PSI JfJ . Although the average density is low, 
the density in the foils is large, so the lower-energy particles 
have a significant probability to stop and decay. At low kinetic 
energies, a positively-charged particle picks up and loses elec- 
trons as it passes through the material [8]. With each passage 
of the beam through a foil, a fraction of the beam exits in the 
neutral charge state (muonium). The muonium atoms are not 
reaccelerated by the electric field and have small probabilities 
to survive. In this way, a large portion of the beam neutralizes 
and the efficiency is low. 

We avoid both problems by using gas as the low-density re- 
tarding medium. Energy is lost and restored over very small dis- 
tances, so that the energy of the [i + fluctuates about r eq . Charge 
exchange processes still proceed in the gas, but muonium atoms 
change back to \i + over short distances, so the (x + are not lost but 
have a reduced effective charge. This decreases the energy gain 
from the electric field. To minimize this decrease, we use he- 
lium gas because in it the effective charge is largest J9[]. 

Elastic scattering off of nuclei places a lower limit on the en- 

i — i 

ergy spread achievable by frictional cooling [2]. The energy 
spread of the cooled beam can be preserved during reacceler- 
ation, yielding higher energy beams with small relative energy 
spreads. 

3. Cooling scheme 

We consider producing a LE-U+ beam using frictional cool- 
ing from a surface muon source such as the uE4 beam at PSI. 

3.1. LE-fi + beam at PSI 

The uE4 beam line delivers the world's highest flux surface 
\x + beam fl. PSI currently produces a LE-(.i + beam by focusing 
the uE4 beam onto a thin foil moderator (about 100 um thick) 
covered with a very thin layer (less than 1 (xm) of a condensed 
van der Waals gas such as argon, neon, or nitrogen cryosolids. 
Very slow (x + exit from the downstream side of the moderator. 
The energy distribution of these muons has a maximum near 
15 eV and a tail extending to higher energies. An electric field 
then reaccelerates the muons yielding a beam with mean energy 
tunable between 0.5 keV to 30 keV with a spread of 400 eV. 
This moderation technique has an efficiency of 10~ 5 to 10 -4 to 
convert a surface muon to a LE-u. + 0]. 

3.2. Frictional cooling scheme 

We searched for a way to efficiently produce a LE-|x + beam 
using frictional cooling. Because frictional cooling only works 
at low energies, we must lower the kinetic energies of the 
muons to below 50 keV. To achieve this, we inject the muons 
into a cooling cell, which is a helium gas cell with an electric 
field in the direction opposite to that of the beam (Fig. [2]). A 
1 T solenoidal magnetic field guides the beam in the cooling 
cell. The gas and the electric field slow the muons down. The 
slowing down is primarily due to the electric field. When they 
turn back, a large fraction of the muons have very low kinetic 
energies and are brought to r eq . This scheme produces a beam 



at the equilibrium energy with a very small energy spread, re- 
gardless of the initial beam energy spread. 
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Figure 2: The block diagram of the frictional cooling scheme. 

To reduce the required electric potential in the cooling cell, 
the beam first passes through a foil for an initial energy loss. 
For a fixed electric field, lowering the mean energy decreases 
the distance that the muons traveled in the cooling cell, thereby 
reducing the loss from decay. 

The cooling cell is ineffective for muons with large transverse 
momentum, because when they turn around they have kinetic 
energies above T' (see Fig.[T). The electric field reaccelerates 
them to approximately their initial energies. Therefore we use 
quadrupoles to reduce the angular divergence of the beam. The 
quadrupoles are designed to maximize the number of muons 
with low transverse momenta. 

Between the quadrupoles and the cooling cell is a weak 
dipole. It has minimal effect on the high-energy input beam, 
letting it pass straight through. After the dipole, the beam en- 
ters the cooling cell. 

A window is needed to separate the cooling cell from the in- 
put beam line. The window must be thin to allow the cooled 
muons get out efficiently. After the muons exit the cooling cell, 
they enter the dipole again. The magnetic field in the y direc- 
tion turns the cooled muons in the x direction. We evaluate the 
momenta and positions of the muons as they exit the dipole. 

4. Simulation 



4.1. Simulation tools 

We have developed a program called CoolSim fioll based on 
Geant4 [11], in which geometries can be implemented easily 
by a macro-command interface. For simulating the frictional 
cooling process, we use the Geant4 low-energy processes Il2ll 
to determine energy loss and spatial displacement. 

We have added new low-energy physics processes to the 
Geant4 framework: hydrogen formation for protons, muonium 
formation for and charge exchange interactions for both 
protons and \i + . The charge exchange process is simulated by 
assigning particles an effective charge, 
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where o~u is the cross section for the changing charge state from 
qi to qj II 1 3fl . The effective charges for \i + in helium and carbon 
are shown in Fig. [3] 
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Table 2: Summary of u,E4 beam properties at PSI 
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Figure 3: The effective charge as a function of kinetic energy for \i* in helium 
(dashed) and carbon (solid). 



We also take the decay process for into account. The en- 
ergy losses calculated by the program were verified against the 
NIST tables. 

We simulated the scheme shown in Fig.[2]in the CoolSim pro- 
gram using uniform electric and magnetic fields without fringe 
fields. The quadrupoles are simulated as a slightly modified 
FODO cell, which is a symmetric ^FODO^F cell. All com- 
ponents have the same aperture. The simulation parameters are 
given in Table Q] 



Table 1 : Parameters of the elements in the cooling scheme 



Aperture 


100 cm 


Foil: 




Thickness 


100 (xm 


Material 


tungsten 


Quadrupoles: 




Gradient 


1 T/m 


Magnet Length 


20 cm 


Drift space 


40 cm 


Dipole: 




Strength 


0.003 T 


Length 


100 cm 


Cooling cell: 




Length 


200 cm 


Electric field strength 


1.8 MV/m 


Magnetic field strength 


1 T 


Gas material 


helium 


Gas density 


0.01 mg/cm 3 


Window: 




Thickness 


20 nm 


Material 


carbon 



4.2. Input beam 

We use the Geant4 General Particle Source to reproduce the 
surface muon beam using the parameters of PSI's |iE4 beam 



Beam energy 
Accepted solid angle 
x/x' (FWHM) 
y/y' (FWHM) 
Ap/p (FWHM) 
u. + rate 



3.7 MeV 
135 msr 

6.5 cm/150 mrad 

2.6 cm/300 mrad 

9.5% 
228 x 10 6 /mAs 




Figure 4: The spatial distribution of the simulated source. 

(Table [2} and assuming gaussian distributions with no correla- 
tions. Figure [4] shows the spatial distribution of the simulated 
source, which agrees well with reference Jj] . The RMS of the 
simulated angular distributions are 0.06 rad in x and 0.12 rad 
in y (Fig. [5]). The kinetic energy distribution of the simulated 
source has a mean value around 3.7 MeV with an RMS of 
0.32 MeV (see Fig.|6]solid curve). 100,000 muons were simu- 
lated. 

4.3. Effect of the foil 

The tungsten foil reduces the mean of the energy distribu- 
tion of the beam from 3.7 MeV to 1.7 MeV (Fig. |6) with 90% 
efficiency. Thicker foils greatly reduce the efficiency, while 
only reducing the mean energy slightly. The angular spreads of 
the beam increase greatly after passing through the foil, which 
makes the quadrupoles more necessary. 

4.4. Effect of the quadrupoles 

The quadrupoles have been designed to rotate the beam phase 
space by 90° (Fig. |7J, enlarging the beam's transverse size and 
increasing the number of muons with low transverse kinetic en- 
ergies. 

Figure [8] shows the longitudinal and transverse kinetic en- 
ergy distributions of the beam directly before the quadrupoles 
and directly before entering the cooling cell. The mean trans- 
verse kinetic energy of the beam is greatly reduced. Figure [8] 
also shows the kinetic energy distribution before entering the 
cooling cell for the subset of the beam that can ultimately be 
cooled and for the subset that cannot be cooled. Muons with 
transverse kinetic energies greater than 50 keV are not cooled. 
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Figure 6: The kinetic energy distribution before (solid) and after (dashed) the 
foil. 



Half of the muons exiting the foil are lost due to the limit of 
the aperture of the quadrupoles. 



4.5. Behavior in the cooling cell 

Figure [9] shows the kinetic energy of representative muons as 
a function of the depth in the cooling cell. The muons enter the 
cooling cell at z = cm with a mean kinetic energy of 1 .7 MeV. 
They are decelerated by the electric field and the gas and drift 
back out towards the entrance of the cell. Four of the muons 
shown in Fig. [9] reach T e „, while two turn back and are reac- 
celerated. The fluctuations in the kinetic energies of the cooled 
muons about T eq are due to large angle scatters off of gas nuclei. 

The depth at which a muon turns around depends on its ki- 
netic energy when it enters the cooling cell. Figure [TOl shows 
the distribution of this depth. The mean depth is 94 cm. A 2 m 
long cooling cell has sufficient length to turn back all the muons 
that can be cooled. The cooled muons have an equilibrium en- 
ergy of 2.2 keV with a spread of 0.26 keV (Fig. QT). 24% of the 
muons from the source are cooled and reach the exit window. 






Figure 8: The transverse (T r ) and longitudinal (7\) kinetic energy distributions 
of the beam directly before the quadrupoles (upper left) and directly before the 
cooling cell (upper right). The lower plots show the distributions directly before 
the cooling cell for the subset of muons that are ultimately cooled (left) and that 
are not cooled (right). 



4.6. Effect of the window 

The window of the cooling cell is 20 nm thick and made of 
carbon. 4% of the cooled muons stop and decay in the window. 
The beam exits the window with a mean energy of 0.9 keV at 
which g e ff = 0.15. When the muons exit the window and enter 
the vacuum, a discrete charge is chosen based on q e ff , which is 
interpreted here as the probability that the muon is not in muo- 
nium atom. 85% of the cooled muons are therefore lost to muo- 
nium formation because the dipole does not turn the muonium 
atoms to the output channel. 

The spatial distribution of the cooled muons exiting the cool- 
ing cell is shown in Fig. [12] The beam has a transverse RMS 
sizes of 24 cm in x and 14 cm in y. 

4.7. Output 

After the dipole, the cooled beam has a mean kinetic en- 
ergy of 0.9 keV and an RMS of 0.32 keV (Fig. [13) ■ Figure [H 
shows a comparison of the energy relative to the mean en- 
ergy for muons at the source and muons exiting the scheme. 
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Figure 9: The kinetic energy as a function of depth in the cooling cell, z = cm 
is the entrance of the cooling cell for representative trajectories. 



Figure 1 1 : The energy distribution of the cooled muons directly before entering 
the window. 





Figure 10: The distribution of maximum depth reached by the cooled muons in Fi § ure 12: The s P atial distribution of the cooled beam exiting the cooling cell 
the cooling cell. ( both ^ + and muomum). 



The peak muon rate is 6.0 x 10 6 /(mAskeV) while the peak 
source rate is 2.9 X 10 s /(mA s keV). The energy spread of the 
beam is greatly reduced. The spatial distribution of muons ex- 
iting the dipole is shown in Fig. [15] The RMS sizes of the 
beam are 20 cm in x and 24 cm in y. 2% of the source muons 
exit the dipole; they would then be collected and reaccelerated. 
From the initial source to the output, the mean time of flight is 
1300 ns (Fig. [HJ. 

4.8. A simplified scheme 

We also simulated a simplified scheme that avoids the com- 
plication of the quadrupoles and greatly reduces the required 
high voltage. 

We limited the design so that the required potential is within 
the available range of common high voltage supplies. The 
quadrupoles are removed but the dipole is kept for extracting 
the cooled muons. The diagram of this scheme is shown in Fig. 
[TTI A foil reduces the energy of the initial beam to maximize 
the number of muons that can be cooled by the cooling cell. 
The thickness of the foil for this scheme is 130 urn. The dipole 



strength is 0.01 T, which turns the cooled muons with a radius 
of 10 cm. The cooling cell is only 30 cm long and has a diam- 
eter of 20 cm. Keeping the electric field at 1.8 MV/m yields 
a potential difference of 540 kV in the cooling cell. AIT 
solenoidal magnetic field guides the beam in the cooling cell. 
The cell is filled with helium gas at a density of 0.01 mg/cm 3 . 
A 20 nm thick carbon window separates the cooling cell from 
the dipole. 

At the output, the cooled beam has a mean kinetic energy of 
0.9 keV and an RMS of 0.32 keV. The overall efficiency in this 
scheme is 0.1%. The peak muon rate is 4.2 x 10 5 /(mAskeV). 

5. Points to investigate 

This paper outlines a frictional cooling scheme for the effi- 
cient production of a LE-(i + beam. There are several issues that 
require further study. 

5.1. Fields 

All fields in the simulations are uniform. The fringe fields 
of a strong solenoid around the cooling cell would change the 
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Figure 13: The energy distribution of the output 
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Figure 15: The transverse distribution of cooled muons exiting the dipole. 
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Figure 14: The energy distribution of the output beam compared to the input 
beam. The distributions are centered on their means. 




t[ns] 

Figure 16: The distribution of time of flight of the cooled muons from source 
to the output. 



muons' trajectories. Bucking coils could be used to reduce the 
fringe fields out of the cooling cell. The dipole can be put fur- 
ther away from the cooling cell to avoid the fringe field of the 
solenoid; or, the fringe field itself can replace the dipole and 
guide the cooled muons out. Realistic scenarios need to be sim- 
ulated. 

An additional radial electric field could focus the beam to 
the axis of the cell, where the radial component of the fringe 
magnetic field is much lower than at the edge of the solenoid. 
How to create such a field should be investigated. 

5.2. Electric fields in gas 

The frictional cooling scheme requires a strong electric field 
in gas. Frictional cooling experiments at the Max-Planck- 
Institute for Physics have used electric fields up to 500 kV/m 
over 10 cm in a gas cell filled with helium at pressures from 
10~ 3 mbar to 1.25 bar without breakdown. Further investiga- 
tions are needed to determine the feasibility of even higher 
fields in helium gas. 

Muons lose a small fraction of their kinetic energy to ion- 
ization in the gas. We estimated an ionization rate of less than 



10 10 /s for a source flux of 4.2 x 10 8 /s. The removal or neu- 
tralization of the resultant charge in the cooling cell will have 
to be investigated. 

5.3. Window 

The window on the cooling cell must be thin and gas tight. 
20 nm thin carbon films can be deposited on nickel grids with 
55% open areaQ. These windows can hold high gas pressures 
over small areas. However, they cut the efficiency of the scheme 
in half. In our case, this would mean 1% overall efficiency. 
The feasibility of large-area thin gas-tight windows needs to 
be investigated. Focusing the beam in the cooling cell would 
be helpful by reducing the area required for the thin window. 
The material of the window should also be optimized to re- 
duce beam neutralization. A layer of noble gas frozen onto 
the window can decrease neutralization of the cooled muons. 
Other technologies such as differential pumping systems could 
be considered to replace the window. 



'Structure Probe Inc.: http://www.2spi.com 
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Figure 17: Block diagram of the simplified scheme. 



6. Conclusion 

Frictional cooling holds promise for producing low energy 
muon beams with small energy spreads. Our simulation shows 
that a frictional cooling scheme can cool a surface muon 
beam to a mean energy of 0.9 keV with an energy spread of 
320 eV. The efficiency for LE-(i + production is greater than 
1%, which is 2 to 3 orders of magnitude higher than the ef- 
ficiency achieved with the moderation technique. If the fric- 
tional cooling scheme is experimentally successful, a surface 
muon rate of 228 x 10 6 /(mA s) would result in a LE-|x + rate of 
2.3 x 10 6 /(mAs). 
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